Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  984 1818 (N42E30 L=174)   960    203        2 ( 941/1756-1027/1880) -32.7 E37-E23 N37-N47
      2  922 1403 (N12E15 L=189)   963    351        0 ( 876/1254- 987/1531) -84.0 E18-E12 N01-N20
      3 1808 1196 (N25W47 L=252)   382     38        0 (1793/1173-1828/1225) -26.6 W45-W49 N23-N27
      4 1125 1169 (N03W02 L=206)   828    248        1 ( 993/1093-1220/1209)   0.0 E05-W08 N01-N05
      5  373  948 (S26E40 L=164)   273     97        0 ( 334/ 924- 417/ 978)  39.8 E43-E37 S28-S23
      6 1023  896 (S14W03 L=208)   640    131        1 ( 959/ 871-1108/ 932)  40.6 W00-W07 S16-S10
      7 1756  886 (N06W50 L=254)   502    126        0 (1715/ 864-1807/ 911)   5.7 W45-W55 N04-N07
      8  770  590 (S39E04 L=200)  1857    509        4 ( 577/ 552-1023/ 625)  17.4 E19-W13 S43-S32
      9 1359  545 (S24W35 L=239)  1213    331        1 (1211/ 516-1511/ 582)  24.2 W25-W45 S28-S19
     10 1489  473 (S24W47 L=251)   348     50        2 (1466/ 445-1516/ 494)  71.6 W45-W48 S25-S22
     11 1198  279 (S45W41 L=246)   311     97        2 (1172/ 261-1225/ 300)  11.3 W38-W44 S46-S43
  

    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