Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  961 1550 (N40W09 L=355)  1664    302        7 ( 851/1501-1092/1611) -36.3 W01-W17 N35-N47
      2 1595 1505 (N21W52 L= 38)   321    100        0 (1564/1471-1634/1524)   0.0 W48-W57 N19-N21
      3 1356 1430 (N23W31 L= 17)   425    142        1 (1337/1372-1395/1479)  50.2 W28-W35 N19-N25
      4 1521 1312 (N12W38 L= 24)   396     76        1 (1492/1278-1552/1346)   9.5 W34-W40 N11-N14
      5  339 1303 (N36E47 L=299)  1077    190        2 ( 251/1256- 422/1336) -13.4 E59-E38 N32-N38
      6  556 1208 (N26E27 L=319)   680    126        3 ( 485/1174- 618/1234) -21.8 E32-E22 N23-N29
      7  482 1026 (N17E35 L=311)   295    108        0 ( 464/ 982- 524/1060)  36.9 E36-E31 N14-N18
      8 1869  995 (S15W61 L= 47)   463     73        0 (1850/ 965-1883/1038)  36.9 W57-W63 S17-S13
      9  304  913 (N12E51 L=295)   433     84        3 ( 277/ 877- 333/ 948) -59.0 E53-E48 N09-N14
     10 1610  846 (S17W33 L= 19)   903    199        4 (1552/ 772-1677/ 934)  35.0 W28-W39 S20-S12
     11  908  655 (S13E14 L=332)   349    139        2 ( 856/ 627- 973/ 702)  14.0 E17-E09 S13-S10
     12 1092  356 (S37E12 L=334)   307     69        1 (1067/ 331-1116/ 387)  31.0 E13-E08 S39-S35
     13  513  326 (S28E59 L=287)   612     77        2 ( 473/ 314- 554/ 345) -14.0 E63-E54 S29-S26
  

    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