Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  389 1489 (N33E50 L=355)   959    272        0 ( 336/1447- 455/1539) -21.8 E57-E42 N30-N36
      2 1564 1514 (N33W44 L= 89)   508     57        1 (1536/1488-1591/1539) -76.0 W41-W45 N31-N35
      3 1742 1515 (N32W61 L=106)   906    151        2 (1706/1490-1774/1539)   0.0 W56-W64 N30-N33
      4  725 1385 (N27E20 L= 25)   265     80        1 ( 685/1370- 770/1400)   9.5 E23-E17 N25-N27
      5 1764 1261 (N16W55 L=101)   428     94        0 (1734/1214-1788/1306) -63.4 W53-W57 N12-N18
      6  258 1219 (N15E57 L=349)  1114    197        1 ( 232/1132- 284/1292) -56.3 E60-E53 N10-N19
      7 1631 1231 (N15W42 L= 88)   467    133        1 (1592/1160-1658/1287) -69.4 W39-W44 N10-N18
      8 1843 1251 (N14W64 L=110)   441     77        0 (1816/1228-1864/1270) -31.0 W61-W66 N12-N15
      9  911 1180 (N14E07 L= 39)   290     81        0 ( 868/1144- 952/1215) -38.7 E09-E04 N11-N15
     10 1314  971 (S00W18 L= 64)   215     74        0 (1274/ 951-1358/ 994)   0.0 W15-W20 S01-N01
     11 1197  570 (S26W11 L= 57)   255     75        0 (1157/ 559-1237/ 585) -11.3 W08-W13 S26-S24
     12  413  547 (S27E47 L=358)   722    150        0 ( 378/ 518- 459/ 569)  24.0 E51-E42 S29-S25
  

    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