Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  334 1539 (N14E64 L=349)  2551    428        4 ( 296/1448- 387/1624) -31.4 E72-E54 N08-N19
      2  814 1487 (N20E24 L= 29)   709    168        1 ( 782/1417- 848/1548) -69.4 E27-E20 N16-N24
      3 1267 1368 (N23W07 L= 60)   635    197        1 (1218/1276-1310/1453) -45.0 W02-W11 N18-N27
      4 1734 1334 (N33W43 L= 96)  2377    268        3 (1613/1216-1863/1437) -15.6 W32-W57 N29-N37
      5  826 1305 (N09E18 L= 35)   405     83        1 ( 800/1272- 860/1343)  78.7 E19-E16 N07-N12
      6  251 1223 (S07E58 L=355)  2060    352        6 ( 164/1131- 350/1301)  36.9 E66-E50 S13-N00
      7 1582 1219 (N22W31 L= 84)   307    135        1 (1545/1173-1616/1243)  63.4 W28-W33 N18-N24
      8  522 1091 (S09E32 L= 21)   260     83        0 ( 479/1073- 558/1108)  21.8 E34-E29 S10-S07
      9 1911  926 (N13W69 L=122)  2381    210        4 (1879/ 869-1943/ 986)  82.9 W65-W75 N08-N16
     10  948  679 (S24W03 L= 56)   362     97        0 ( 896/ 665-1006/ 692)  26.6 W00-W06 S26-S22
     11 1467  293 (S35W56 L=109)   704     87        1 (1432/ 274-1501/ 311)   8.1 W52-W59 S37-S33
     12 1381  274 (S39W51 L=104)   401     39        0 (1363/ 257-1399/ 297) -71.6 W49-W53 S40-S37
  

    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