Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  795 1676 (N51W05 L=264)   435     84        1 ( 770/1639- 814/1717)  56.3 W03-W09 N47-N53
      2 1462 1686 (N29W54 L=312)   363     46        0 (1442/1662-1485/1709)   8.1 W50-W57 N28-N30
      3  411 1260 (N33E35 L=223)   372    156        0 ( 352/1237- 495/1289) -21.0 E41-E28 N31-N36
      4  307 1194 (N31E46 L=213)   721    122        1 ( 254/1156- 359/1230)   5.2 E51-E40 N29-N32
      5  639 1085 (N17E21 L=238)   648    119        1 ( 613/1032- 664/1150)  71.6 E23-E18 N14-N20
      6  462 1087 (N21E33 L=225)   448    128        1 ( 404/1049- 524/1114) -32.0 E37-E29 N17-N24
      7 1701  710 (S33W38 L=296)   391     66        0 (1672/ 690-1735/ 736)   9.5 W34-W41 S33-S31
      8 1412  651 (S27W14 L=272)   298     90        1 (1369/ 621-1438/ 683)  11.3 W10-W15 S28-S25
      9  528  427 (S18E50 L=208)   362     65        1 ( 507/ 391- 547/ 449)  26.6 E53-E47 S20-S16
     10 1276  197 (S60E13 L=246)   601     75        1 (1254/ 176-1311/ 233) -82.9 E17-E09 S63-S54
     11 1193  172 (S59E23 L=236)   947     80        1 (1163/ 140-1225/ 203)  82.9 E29-E18 S62-S53
     12 1096  153 (S56E33 L=225)   590     64        0 (1070/ 137-1121/ 176)  -9.5 E37-E29 S58-S53
  

    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