Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  417 1606 (N15E61 L= 17)   461     74        0 ( 400/1577- 445/1641) -24.0 E65-E56 N13-N17
      2  287 1514 (N07E70 L=  8)  1283    115        1 ( 249/1462- 323/1556)  63.4 E74-E67 N03-N09
      3 1336 1501 (N33W06 L= 84)   346     70        0 (1302/1471-1371/1523) -11.3 W03-W08 N31-N34
      4 1897 1014 (N21W65 L=144)  2385    284        0 (1832/ 907-1949/1125)  42.0 W58-W80 N12-N30
      5  656 1011 (S14E21 L= 57)   241     79        0 ( 610/ 998- 694/1025)   9.5 E24-E18 S15-S13
      6  318  817 (S34E46 L= 33)   362    119        0 ( 278/ 784- 374/ 836)  20.0 E50-E39 S37-S32
      7  562  698 (S35E22 L= 57)  1832    365        5 ( 459/ 588- 689/ 830)  59.4 E30-E14 S45-S23
      8 1608  591 (S12W51 L=129)  1189    289        1 (1496/ 482-1667/ 689)  82.9 W48-W55 S20-S04
      9 1281  482 (S28W34 L=112)   942    238        2 (1176/ 425-1370/ 535)  59.7 W29-W38 S34-S22
     10 1132  414 (S37W28 L=106)   454     89        2 (1102/ 385-1166/ 451) -15.9 W24-W31 S38-S34
  

    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