Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  806 1619 (N22E29 L=240)   585    105        0 ( 755/1591- 871/1641)  45.0 E32-E25 N19-N24
      2  349 1584 (N11E64 L=206)   648     52        2 ( 325/1558- 379/1607)  45.0 E68-E60 N09-N12
      3  354 1404 (N00E53 L=217)  1061    286        3 ( 276/1312- 439/1519)  65.2 E57-E49 S05-N08
      4 1671 1430 (N35W33 L=303)  1467    265        1 (1556/1351-1775/1495)  -5.2 W22-W44 N33-N37
      5  319 1235 (S10E50 L=219)   584    110        1 ( 280/1206- 356/1268) -14.0 E54-E46 S11-S08
      6  555 1192 (S08E31 L=238)   536    187        0 ( 532/1108- 597/1263) -77.5 E33-E28 S12-S03
      7  917 1198 (S00E10 L=260)   214     81        0 ( 901/1160- 936/1240) -33.7 E11-E08 S01-N01
      8 1910 1095 (N23W62 L=332)   994    122        1 (1883/1040-1928/1152)  76.0 W58-W65 N18-N26
      9 1016 1132 (S02E03 L=267)   210     30        1 ( 993/1111-1039/1150)  45.0 E04-E01 S02-N00
     10 1841 1067 (N18W53 L=323)   373     53        0 (1821/1048-1864/1095)  90.0 W50-W55 N16-N19
     11 1171  720 (S20W16 L=286)   963    288        2 (1129/ 645-1217/ 862) -63.4 W13-W21 S23-S11
     12  654  517 (S47E12 L=258)  1232    257        2 ( 544/ 468- 766/ 576)   0.0 E22-E01 S48-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