Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  783 1650 (N30E34 L= 97)   520    104        2 ( 726/1610- 811/1681)  -7.1 E38-E30 N27-N32
      2  515 1522 (N16E47 L= 84)   890    184        3 ( 454/1482- 571/1559)  18.4 E51-E41 N12-N18
      3  787 1449 (N17E25 L=106)   558    184        1 ( 726/1406- 860/1471)  32.0 E28-E20 N14-N20
      4 1382 1383 (N27W14 L=145)   500     84        1 (1352/1337-1400/1425) -71.6 W12-W15 N24-N30
      5  888 1390 (N16E17 L=114)   259     61        1 ( 859/1372- 915/1417)  14.0 E18-E14 N14-N17
      6  551 1161 (S04E32 L= 99)   277     91        0 ( 523/1121- 585/1190)  68.2 E33-E29 S07-S02
      7  524  861 (S23E28 L=103)  1145    376        4 ( 409/ 765- 625/1007)  55.6 E37-E22 S31-S12
      8 1361  528 (S24W35 L=166)   661     96        3 (1325/ 498-1403/ 570)  63.4 W32-W38 S25-S20
      9  706  449 (S46E08 L=123)  1301    257        4 ( 588/ 426- 825/ 475)  16.7 E19-W01 S49-S41
     10 1088  181 (S57W44 L=175)  1107    103        2 (1051/ 142-1121/ 229) -39.8 W38-W54 S61-S51
  

    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