Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1061 1562 (N39W07 L=305)   230     47        0 (1035/1547-1088/1577)   0.0 W04-W08 N38-N40
      2 1170 1165 (N12W10 L=308)   301    139        0 (1112/1135-1228/1183) -15.9 W06-W13 N10-N13
      3  188 1109 (N12E64 L=234)   502     85        0 ( 167/1071- 211/1150) -33.7 E67-E60 N10-N14
      4  854  890 (S02E11 L=287)  1564    486        1 ( 781/ 683- 964/1101) -85.4 E14-E04 S15-N10
      5  522 1047 (N08E32 L=266)   217     75        0 ( 511/1003- 536/1079)  78.7 E33-E31 N05-N10
      6 1671  968 (S04W43 L=341)   713    199        0 (1593/ 923-1738/1042)  18.4 W36-W49 S05-N00
      7 1688  830 (S12W44 L=342)   405    115        0 (1668/ 774-1723/ 871)  51.3 W42-W47 S15-S10
      8  400  663 (S15E46 L=252)   477     95        2 ( 374/ 620- 427/ 707)  80.5 E48-E44 S18-S12
      9 1457  588 (S26W27 L=325)   353     96        0 (1422/ 556-1499/ 624)  31.0 W24-W30 S27-S23
     10 1242  346 (S43W13 L=310)  1176    155        5 (1194/ 288-1287/ 412) -80.5 W08-W16 S48-S36
     11  897  290 (S45E17 L=280)   325     51        1 ( 875/ 267- 913/ 319)  51.3 E19-E15 S47-S42
     12  702  252 (S48E39 L=259)  2452    286        4 ( 570/ 213- 830/ 293)   4.8 E64-E24 S52-S44
  

    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