Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  233 1288 (S10E62 L=271)   688     96        1 ( 206/1238- 265/1336)  74.1 E63-E59 S13-S06
      2  187 1095 (S22E60 L=272)   805    163        0 ( 149/1042- 224/1148)  54.5 E63-E57 S25-S18
      3  569 1100 (S14E29 L=304)   208     46        0 ( 557/1075- 586/1133) -63.4 E29-E27 S15-S11
      4 1603 1050 (N12W33 L=  6)   258     73        1 (1590/1010-1616/1090) -53.1 W31-W34 N10-N14
      5  519  946 (S24E28 L=304)   557    106        0 ( 470/ 915- 576/ 984)  -6.3 E32-E23 S24-S22
      6  224  856 (S36E53 L=280)   460    111        0 ( 199/ 823- 257/ 880)  38.7 E55-E49 S38-S34
      7 1739  724 (S00W55 L= 27)   414     62        1 (1717/ 695-1761/ 752) -56.3 W53-W56 S02-N01
      8 1153  476 (S33W27 L=359)   431    100        1 (1123/ 442-1188/ 518) -14.0 W22-W30 S33-S30
      9 1617  461 (S18W59 L= 31)   383     95        0 (1582/ 445-1665/ 486)  45.0 W55-W61 S19-S14
     10 1355  388 (S31W44 L= 17)   964    185        3 (1280/ 359-1426/ 427)  15.3 W38-W49 S32-S28
  

    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