Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1159 1685 (N48W08 L=251)   273     49        0 (1140/1661-1187/1701)  36.9 W06-W10 N45-N49
      2 1295 1646 (N45W20 L=263)   440     88        0 (1256/1621-1340/1666)  20.6 W16-W24 N43-N46
      3 1566 1603 (N43W46 L=289)  2000    273        3 (1432/1563-1680/1633)  -3.9 W31-W62 N39-N44
      4  518 1380 (N22E37 L=206)   278     46        1 ( 493/1360- 552/1403)  21.8 E39-E34 N20-N23
      5 1561 1127 (N11W35 L=278)  2211    415        4 (1484/ 943-1615/1294) -76.0 W29-W38 N00-N20
      6  803 1202 (N12E15 L=228)   290    115        1 ( 772/1159- 843/1249) -50.2 E17-E11 N09-N15
      7 1100 1215 (N15W04 L=247)   222     97        0 (1059/1198-1147/1242)  31.0 W01-W06 N13-N16
      8  252 1079 (N01E55 L=188)  1491    217        3 ( 187/ 996- 291/1179)  55.0 E62-E51 S03-N07
      9 1023 1121 (N08E01 L=242)   337    109        0 (1001/1078-1050/1178) -80.5 E01-W01 N05-N11
     10  846  993 (S00E11 L=232)   904    255        2 ( 787/ 882- 892/1085)  85.2 E14-E07 S07-N05
     11 1088  718 (S16W05 L=248)  2037    622        4 ( 992/ 485-1211/ 971)  71.0 W00-W11 S32-N00
     12  442  866 (S10E38 L=205)   880    235        0 ( 400/ 779- 515/ 965) -54.0 E41-E32 S15-S04
     13 1643  784 (S10W43 L=286)   342     79        0 (1625/ 738-1667/ 822)  80.5 W42-W44 S13-S07
     14 1733  705 (S15W53 L=296)   996    120        2 (1698/ 654-1785/ 754)  35.5 W50-W58 S18-S11
     15 1655  649 (S19W47 L=290)   503    124        0 (1614/ 619-1706/ 699) -21.8 W42-W52 S20-S15
     16 1344  277 (S48W36 L=279)   405     41        1 (1317/ 263-1370/ 295)   0.0 W33-W40 S49-S46
     17 1037  260 (S51W06 L=249)   461     66        1 ( 997/ 242-1073/ 277)  15.9 W02-W09 S53-S49
     18  790  245 (S55E20 L=223)  1779    186        3 ( 709/ 222- 884/ 276)   5.7 E30-E09 S58-S51
  

    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