Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  931 1547 (N28E09 L=224)   263     70        1 ( 911/1509- 950/1577) -53.1 E10-E07 N25-N29
      2  810 1488 (N23E16 L=217)   545    138        0 ( 723/1458- 867/1504) -11.3 E22-E12 N21-N24
      3  716 1412 (N18E22 L=212)   237     50        0 ( 684/1391- 746/1432) -26.6 E23-E19 N16-N19
      4  605 1365 (N15E28 L=205)   312     62        1 ( 591/1325- 622/1401)  90.0 E29-E27 N12-N16
      5 1901 1219 (N14W65 L=298)  2005    355        1 (1833/1132-1970/1316) -22.8 W57-W77 N09-N19
      6 1215 1263 (N11W10 L=244)   420     93        0 (1175/1221-1253/1309) -45.0 W07-W12 N08-N13
      7 1591 1214 (N11W35 L=269)   491    166        1 (1526/1195-1661/1251)   0.0 W30-W40 N09-N12
      8  411 1181 (N03E40 L=193)   230     41        0 ( 394/1159- 437/1197)  18.4 E41-E38 N01-N03
      9  459 1043 (S05E36 L=197)   792    143        2 ( 436/ 960- 497/1113) -83.7 E37-E33 S10-S01
     10 1869  981 (S01W62 L=295)   418     60        0 (1852/ 949-1890/1006)  71.6 W59-W64 S02-N01
     11  157  866 (S16E65 L=169)  2466    371        6 ( 125/ 767- 196/ 968)  67.4 E69-E60 S21-S09
     12 1599  909 (S08W37 L=271)   244     22        1 (1575/ 888-1618/ 927)  18.4 W35-W38 S08-S06
     13  842  886 (S13E10 L=223)   252     58        0 ( 817/ 865- 879/ 906)   0.0 E12-E08 S14-S12
     14 1153  760 (S20W09 L=243)  1242    293        0 (1043/ 657-1319/ 843)  39.1 W03-W19 S26-S13
     15 1795  545 (S27W64 L=297)  1267    168        0 (1738/ 507-1849/ 605)  37.9 W59-W68 S30-S22
     16  960  548 (S34E02 L=231)   217     48        0 ( 941/ 530- 985/ 570) -18.4 E03-W00 S35-S32
     17  883  498 (S38E07 L=226)   255     81        0 ( 842/ 476- 920/ 528)  26.6 E10-E04 S39-S35
  

    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