Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  703 1687 (N38E28 L= 81)  3904    736       12 ( 503/1587- 942/1785) -17.7 E52-E08 N30-N46
      2  754 1417 (N19E19 L= 90)   255    139        0 ( 708/1363- 803/1461) -40.6 E22-E15 N15-N21
      3  606 1340 (N13E28 L= 82)   361     84        0 ( 592/1301- 622/1386) -80.5 E29-E26 N10-N16
      4 1318 1295 (N13W17 L=126)   272     82        1 (1280/1270-1360/1332) -21.8 W14-W19 N11-N15
      5  966 1021 (S05E03 L=106)   232     61        1 ( 938/ 987-1002/1052)  45.0 E04-E01 S06-S02
      6  794  695 (S25E14 L= 96)  1000    290        1 ( 677/ 600- 929/ 796)  37.4 E22-E04 S31-S18
      7 1570  726 (S19W38 L=147)   264     35        0 (1545/ 703-1596/ 751)  45.0 W36-W39 S20-S17
      8  246  679 (S27E57 L= 52)   674    141        0 ( 226/ 632- 271/ 728)  45.0 E60-E53 S29-S23
      9 1427  668 (S24W29 L=138)   335    109        1 (1362/ 653-1475/ 688)  15.9 W24-W32 S24-S22
     10 1692  664 (S22W49 L=159)   332     55        1 (1666/ 647-1717/ 684) -21.8 W46-W51 S22-S20
     11 1103  549 (S33W08 L=117)   510    126        1 (1025/ 532-1164/ 567)   0.0 W02-W12 S34-S32
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y