Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  322 1537 (N13E65 L=350)  2346    394        7 ( 288/1446- 373/1618) -30.5 E73-E56 N08-N18
      2  791 1492 (N20E26 L= 29)   754    187        4 ( 761/1431- 833/1559) -36.9 E29-E21 N16-N23
      3 1238 1378 (N23W05 L= 60)   524    183        2 (1195/1284-1281/1456) -45.0 W00-W09 N18-N27
      4 1712 1343 (N34W41 L= 96)  1900    353        4 (1600/1231-1831/1442) -20.0 W30-W53 N28-N37
      5  793 1314 (N09E20 L= 35)   469     65        3 ( 763/1277- 824/1348) -56.3 E22-E18 N07-N11
      6  240 1238 (S06E60 L=355)  1660    307        5 ( 172/1155- 328/1303)  45.0 E66-E52 S11-N00
      7  406 1110 (S10E41 L= 14)  1279    327        3 ( 299/1078- 539/1147)   6.0 E50-E31 S12-S07
      8 1903  931 (N13W68 L=122)  2406    244        5 (1864/ 870-1937/1000)  77.5 W63-W72 N08-N17
      9  922  692 (S24W01 L= 56)   310     68        2 ( 881/ 675- 964/ 709)  11.3 E01-W04 S25-S22
     10 1487  448 (S25W47 L=102)   320     77        1 (1454/ 418-1504/ 484)  90.0 W45-W49 S27-S22
     11 1454  299 (S35W55 L=110)   777     90        2 (1417/ 279-1485/ 325)  78.7 W51-W58 S37-S32
     12 1372  280 (S39W50 L=105)   407     55        1 (1353/ 259-1391/ 306) -45.0 W48-W52 S40-S37
  

    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