Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  993 1613 (N33E05 L=293)   572     84        1 ( 968/1558-1021/1659) -82.9 E06-E02 N28-N36
      2  581 1623 (N33E36 L=262)   366     65        0 ( 545/1609- 621/1640)   9.5 E38-E32 N31-N34
      3 1433 1478 (N26W26 L=323)   712    161        0 (1359/1441-1521/1521) -14.0 W20-W32 N23-N28
      4 1626 1302 (N16W39 L=336)   736    224        1 (1527/1234-1708/1345) -17.1 W31-W44 N11-N17
      5 1892 1212 (N12W64 L=  2)  1981    223        3 (1863/1125-1945/1318) -65.6 W61-W72 N07-N19
      6 1760 1249 (N13W50 L=347)   345     52        0 (1742/1215-1781/1278) -63.4 W48-W51 N11-N15
      7  176 1248 (N08E64 L=234)   717     72        0 ( 133/1224- 216/1271)  -8.7 E71-E58 N06-N10
      8 1588 1187 (N08W35 L=333)   299     93        0 (1540/1167-1636/1217)  23.2 W31-W38 N06-N10
      9 1617 1035 (S01W37 L=335)   582    147        0 (1591/ 963-1653/1107) -54.5 W35-W40 S04-N03
     10 1179  851 (S14W10 L=308)   255    103        0 (1161/ 801-1195/ 904)  90.0 W09-W10 S16-S10
     11 1858  819 (S11W63 L=  0)  2047    162        8 (1800/ 747-1911/ 871)  32.5 W56-W69 S15-S07
     12 1016  641 (S27W01 L=299)  1643    451        1 ( 917/ 517-1147/ 751)  34.7 E05-W08 S35-S20
     13 1668  589 (S27W51 L=348)  1090    154        1 (1587/ 560-1735/ 619)  13.0 W42-W56 S28-S24
     14  268  574 (S33E63 L=235)   745     94        1 ( 230/ 542- 313/ 597)   3.8 E71-E56 S34-S31
     15 1484  372 (S43W45 L=342)  2855    369        4 (1312/ 341-1665/ 414)   7.5 W27-W65 S45-S38
     16  671  339 (S50E30 L=267)  1990    294        3 ( 506/ 311- 801/ 384)   3.4 E50-E16 S52-S46
  

    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