Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1208 1710 (N40E06 L= 21)   526     84        1 (1169/1684-1242/1738)  45.0 E08-E03 N37-N42
      2 1298 1634 (N37W03 L= 30)   323     31        0 (1275/1610-1323/1659) -36.9 W00-W04 N35-N38
      3  428 1551 (N13E53 L=334)  1154    237        1 ( 402/1461- 448/1647) -65.2 E57-E49 N06-N19
      4  968 1411 (N14E12 L= 15)   657    179        0 ( 939/1324-1019/1491) -48.4 E14-E06 N09-N18
      5 1357 1460 (N26W10 L= 37)   225     41        0 (1336/1433-1378/1482) -26.6 W08-W12 N25-N27
      6  312 1262 (S07E52 L=335)   340     65        0 ( 289/1229- 337/1299)  63.4 E53-E50 S08-S04
      7 1785 1182 (N21W46 L= 73)   798    200        0 (1733/1099-1838/1265) -26.6 W40-W52 N17-N25
      8  533 1106 (S12E31 L=356)   815    178        1 ( 474/1027- 590/1193)  65.6 E34-E28 S17-S06
      9  144 1116 (S17E66 L=321)   605    113        2 ( 127/1064- 158/1166) -45.0 E70-E63 S19-S13
     10 1140  969 (S08W08 L= 35)  1101    293        2 (1039/ 833-1253/1101)  72.5 W05-W11 S17-N02
     11  846  803 (S24E06 L= 21)   662    141        1 ( 773/ 782- 921/ 828)  24.0 E10-E01 S25-S20
     12 1667  714 (S08W49 L= 76)  1149    193        1 (1577/ 656-1731/ 785)  57.5 W43-W52 S13-S02
  

    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