Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  940 1616 (N37E16 L=354)   407    139        0 ( 896/1571- 993/1667) -24.4 E21-E10 N34-N40
      2 1279 1544 (N37W11 L= 21)   300     51        1 (1253/1528-1311/1569) -14.0 W08-W13 N35-N39
      3 1743 1444 (N37W56 L= 66)   457     66        0 (1704/1430-1775/1459)   5.2 W50-W61 N36-N38
      4  514 1374 (N15E39 L=331)   586    194        0 ( 484/1276- 536/1456) -48.0 E44-E35 N09-N19
      5  242 1372 (N11E66 L=304)   883    101        2 ( 220/1333- 274/1423) -50.2 E69-E62 N08-N14
      6  760 1222 (N09E19 L=351)   876    231        3 ( 715/1111- 813/1313)  85.9 E20-E16 N01-N15
      7  440 1118 (S02E39 L=331)   427     68        0 ( 410/1086- 477/1157) -21.8 E42-E36 S03-N00
      8  981  663 (S22W02 L= 13)  1981    650        1 ( 729/ 564-1247/ 777)  32.8 E14-W17 S32-S11
      9  335  600 (S36E51 L=319)   970    212        0 ( 293/ 529- 395/ 692) -59.0 E54-E46 S41-S30
     10  491  628 (S32E34 L=337)   378     81        0 ( 450/ 617- 537/ 645)   0.0 E37-E29 S32-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