Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  782 1660 (N48E00 L=220)   285     31        1 ( 757/1636- 799/1679)   0.0 E02-W01 N45-N49
      2 1586 1613 (N23W56 L=276)  1020    158        4 (1532/1559-1644/1666) -54.5 W52-W59 N19-N26
      3  588 1624 (N51E19 L=202)   404     71        0 ( 553/1599- 612/1655)  31.0 E22-E16 N49-N53
      4 1229 1597 (N31W29 L=249)  1971    412        7 (1051/1542-1420/1633) -27.6 W16-W39 N23-N35
      5 1680 1529 (N16W59 L=279)   954    177        4 (1617/1464-1728/1588)  15.5 W49-W67 N13-N18
      6 1554 1518 (N19W47 L=268)   286     85        1 (1511/1499-1588/1538) -45.0 W44-W49 N16-N20
      7 1662 1328 (N06W47 L=268)   930    191        2 (1580/1299-1717/1361) -33.7 W41-W52 N02-N08
      8  593 1162 (N20E23 L=198)  1276    398        7 ( 405/1102- 751/1229) -27.6 E36-E13 N13-N27
      9  127 1103 (N25E68 L=152)   729     69        1 (  97/1074- 159/1128) -21.8 E75-E63 N23-N27
     10  321  951 (N13E47 L=174)   499    110        4 ( 287/ 919- 361/ 994) -33.7 E49-E42 N11-N16
     11 1428  931 (S11W21 L=242)   333     84        2 (1375/ 903-1458/ 957) -21.8 W18-W23 S12-S09
     12 1399  824 (S16W17 L=238)   334    114        1 (1361/ 801-1441/ 845) -11.3 W15-W20 S18-S15
     13  867  697 (S12E16 L=204)   530    160        0 ( 838/ 640- 907/ 763)  76.0 E17-E13 S16-S08
     14 1875  739 (S35W61 L=282)  1053    120        0 (1816/ 718-1934/ 756)  -7.1 W52-W76 S36-S33
     15 1479  633 (S30W20 L=241)   707    190        0 (1390/ 578-1531/ 702)  15.3 W13-W24 S32-S26
     16  840  436 (S27E26 L=194)  1915    462        3 ( 667/ 356-1044/ 551)   5.5 E41-E10 S29-S24
     17 1222  339 (S44E05 L=216)  1764    364        1 (1055/ 283-1361/ 436)   2.2 E19-W08 S47-S40
     18  400  392 (S21E64 L=156)   610     33        0 ( 380/ 357- 419/ 418)  26.6 E69-E60 S23-S19
  

    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