Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  984 1697 (N41E21 L= 89)   374     48        1 ( 954/1676-1019/1712)  26.6 E22-E18 N39-N42
      2  377 1641 (N24E71 L= 39)  3408    155        5 ( 323/1583- 422/1700) -71.6 E87-E64 N19-N28
      3  983 1400 (N20E10 L=100)   473    108        1 ( 941/1352-1016/1455) -50.2 E13-E07 N17-N23
      4  385 1281 (N02E48 L= 62)   589    140        2 ( 338/1224- 439/1323)  50.2 E50-E43 S02-N05
      5 1028 1300 (N15E05 L=105)   325     67        0 ( 994/1278-1061/1322)  36.9 E07-E03 N13-N16
      6  650 1094 (S04E24 L= 86)   766    156        2 ( 587/1052- 713/1122)  40.6 E27-E20 S08-S01
      7 1862  977 (N13W63 L=173)   473     77        0 (1847/ 937-1874/1015)  90.0 W61-W64 N10-N15
      8 1069  837 (S11W06 L=116)   865    331        1 ( 944/ 730-1180/ 926)  60.6 W00-W11 S20-S04
      9  327  888 (S23E46 L= 64)   323     85        1 ( 292/ 861- 364/ 905)  29.7 E49-E42 S25-S21
     10 1524  786 (S05W36 L=146)   264     84        0 (1494/ 760-1546/ 820) -71.6 W34-W38 S06-S03
     11 1414  687 (S13W32 L=142)   465    114        0 (1384/ 635-1458/ 729)  71.6 W30-W33 S16-S10
     12 1102  638 (S23W13 L=123)   209     86        0 (1082/ 601-1122/ 677)  90.0 W12-W14 S25-S20
  

    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