Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  972 1822 (N42E31 L=175)  1003    225        1 ( 931/1762-1015/1890) -31.0 E39-E24 N37-N47
      2  887 1369 (N09E16 L=189)   968    391        1 ( 772/1225- 963/1531)  87.3 E19-E13 S01-N20
      3 1046 1203 (N03E03 L=203)   371    130        2 ( 971/1182-1117/1217)  20.6 E07-W01 N01-N05
      4  748  920 (S20E13 L=192)   247     93        0 ( 728/ 894- 792/ 958) -36.9 E15-E10 S20-S17
      5 1009  908 (S14W02 L=208)   549    130        1 ( 941/ 870-1078/ 945)  39.8 E01-W05 S17-S10
      6 1736  897 (N06W48 L=253)   302    110        1 (1703/ 877-1780/ 918)   7.1 W44-W52 N04-N06
      7  726  599 (S40E08 L=198)  1447    324        3 ( 572/ 575- 890/ 633)  14.0 E20-W04 S43-S34
      8 1338  554 (S24W33 L=239)  1060    331        0 (1184/ 527-1490/ 591)  29.1 W23-W42 S29-S19
  

    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