Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  862 1657 (N29E29 L=233)  1426    319        5 ( 773/1532- 927/1732) -24.2 E39-E19 N22-N33
      2  662 1627 (N22E41 L=221)  1364    242        5 ( 582/1540- 735/1711)  85.9 E43-E37 N15-N29
      3 1354 1395 (N26W10 L=272)  1863    484        1 (1243/1176-1408/1560) -49.6 W00-W17 N13-N35
      4 1259 1177 (N10W09 L=271)   201     69        0 (1230/1153-1288/1202) -14.0 W07-W11 N09-N11
      5 1657 1025 (N13W39 L=301)   301     85        0 (1642/ 985-1667/1078) -68.2 W37-W39 N10-N15
      6  293 1025 (S21E49 L=213)   713    171        1 ( 237/ 966- 363/1070)  42.0 E54-E42 S25-S16
      7  682  882 (S20E17 L=245)   215     90        1 ( 656/ 852- 713/ 930) -45.0 E18-E14 S21-S17
      8  309  905 (S28E46 L=216)   494    107        1 ( 259/ 894- 370/ 920)  15.3 E51-E40 S28-S25
      9 1161  864 (S09W12 L=274)   345    137        0 (1095/ 831-1220/ 892)   7.1 W07-W15 S10-S08
     10 1781  773 (N02W57 L=319)   686    132        3 (1731/ 730-1824/ 830)  60.3 W53-W60 S01-N06
     11  644  437 (S51E11 L=251)  1644    299        3 ( 524/ 363- 760/ 525) -11.3 E23-W02 S54-S47
  

    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