Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  765 1543 (N30E18 L=  3)   692    185        2 ( 707/1491- 833/1598) -38.7 E22-E12 N25-N33
      2  461 1059 (S00E36 L=345)   483    153        0 ( 394/1026- 516/1095) -12.5 E40-E31 S02-N02
      3  432  881 (S11E39 L=342)   922    211        2 ( 407/ 774- 456/ 996)  74.1 E41-E36 S17-S03
      4  800  857 (S13E14 L=  6)   188     92        0 ( 755/ 840- 847/ 872)   0.0 E16-E10 S13-S12
      5 1396  824 (S15W23 L= 44)   276    116        0 (1333/ 805-1445/ 851)  15.9 W19-W26 S16-S13
      6  213  790 (S15E60 L=321)  1224    208        1 ( 180/ 728- 242/ 850) -76.0 E63-E56 S19-S11
      7  784  662 (S25E16 L=  4)   377     47        0 ( 764/ 626- 807/ 698)  68.2 E17-E14 S27-S22
      8 1103  643 (S27W05 L= 25)   271    104        0 (1076/ 598-1116/ 697)  80.5 W03-W05 S29-S23
      9  502  579 (S30E38 L=342)   718    203        2 ( 456/ 529- 570/ 645)  30.3 E43-E31 S32-S25
     10  257  542 (S31E62 L=318)   586    138        3 ( 236/ 515- 281/ 568)  21.8 E65-E59 S32-S29
  

    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