Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1551 1617 (N32W53 L= 33)   440    141        0 (1515/1590-1577/1650)   7.1 W48-W56 N30-N33
      2 1283 1494 (N30W27 L=  7)   812    193        3 (1214/1447-1370/1534)  -4.4 W21-W34 N28-N32
      3  371 1478 (N44E47 L=294)   359     96        0 ( 341/1455- 410/1492)   9.5 E49-E42 N43-N45
      4 1761 1428 (N16W63 L= 43)  1242    198        2 (1703/1401-1815/1454) -17.1 W56-W69 N12-N18
      5  685 1131 (N18E20 L=321)   661    233        0 ( 574/1075- 780/1195) -39.8 E26-E14 N13-N23
      6 1369  962 (S03W20 L=  0)  1918    455        5 (1314/ 766-1409/1157)  74.7 W13-W22 S14-N09
      7 1603  989 (S06W36 L= 17)   357     98        2 (1575/ 945-1637/1019) -90.0 W34-W38 S08-S03
      8  143  936 (N10E71 L=270)   642    224        1 ( 107/ 898- 167/ 989)  23.2 E79-E65 N07-N13
      9 1545  721 (S21W29 L= 10)  1613    390        0 (1497/ 589-1576/ 911)  84.6 W26-W31 S30-S08
     10 1121  373 (S38E06 L=335)  1522    348        3 (1005/ 281-1191/ 480)  35.8 E17-W01 S43-S30
  

    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