Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  862 1762 (N45E32 L=209)   615     45        1 ( 839/1727- 892/1792) -50.2 E34-E28 N41-N47
      2  278 1518 (N17E71 L=170)   888     58        1 ( 250/1497- 313/1542)  -3.6 E80-E64 N15-N18
      3 1747 1489 (N42W55 L=296)   572     79        0 (1706/1466-1791/1507)   0.0 W49-W62 N41-N43
      4 1256 1327 (N21W09 L=250)   234    109        1 (1213/1286-1286/1363) -31.0 W06-W11 N19-N23
      5 1082 1243 (N13E00 L=241)   885    276        1 ( 991/1196-1166/1283)  24.0 E05-W04 N09-N16
      6  811 1157 (N04E15 L=226)   445    184        0 ( 796/1077- 825/1250) -78.7 E16-E13 S01-N09
      7 1079 1173 (N09W01 L=242)   764    386        1 ( 980/1122-1239/1227)   7.6 E05-W10 N05-N13
      8 1320 1078 (N07W17 L=258)   221    102        0 (1279/1060-1365/1095)  11.3 W14-W19 N05-N08
      9  183  902 (S22E63 L=178)  1332    204        1 ( 119/ 844- 240/ 975)  26.6 E77-E54 S26-S16
     10 1064  839 (S11W05 L=246)   847    312        2 (1001/ 711-1134/ 964)  76.8 W02-W07 S19-S02
     11 1497  656 (S15W37 L=278)   395     95        1 (1450/ 623-1536/ 685) -15.9 W33-W40 S16-S13
     12  531  356 (S56E34 L=206)  1533    169        2 ( 469/ 312- 604/ 402) -18.4 E42-E25 S59-S51
  

    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