Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1472 1640 (N48W13 L=357)   710    139        1 (1414/1614-1516/1670)  41.2 W08-W17 N43-N52
      2  389 1572 (N11E60 L=283)   570    122        1 ( 361/1547- 419/1602) -23.2 E63-E56 N10-N13
      3  794 1568 (N20E30 L=314)   265     60        1 ( 766/1540- 816/1596)  90.0 E30-E28 N18-N22
      4 1195 1256 (N13W03 L=347)  1399    408        3 (1096/1047-1308/1436) -40.4 E06-W14 N03-N20
      5 1713 1289 (N31W38 L= 22)   302     64        0 (1687/1259-1737/1309) -26.6 W35-W40 N29-N32
      6  255 1217 (S12E57 L=287)   431     95        1 ( 222/1202- 287/1235)  15.9 E60-E53 S13-S11
      7  191 1183 (S15E63 L=281)   644    111        1 ( 162/1157- 222/1203)   0.0 E66-E58 S16-S13
      8 1822 1047 (N20W52 L= 36)   394     69        0 (1804/1013-1841/1083) -31.0 W49-W54 N18-N21
      9 1746  904 (N09W49 L= 33)   400     79        0 (1714/ 868-1768/ 943)  68.2 W47-W50 N06-N11
     10 1268  733 (S14W22 L=  6)   478    178        0 (1181/ 695-1351/ 786)  52.1 W18-W25 S18-S08
     11 1729  579 (S08W62 L= 46)   990    162        1 (1665/ 539-1775/ 640)  60.9 W57-W65 S12-S03
     12 1477  420 (S25W50 L= 34)  1284    222        2 (1381/ 391-1568/ 453)  23.2 W42-W56 S27-S20
     13 1228  261 (S43W45 L= 29)   354     40        0 (1203/ 243-1249/ 277) -11.3 W42-W47 S44-S41
  

    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