Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  842 1511 (N41E00 L=254)   288     63        0 ( 805/1497- 876/1527)   0.0 E03-W02 N39-N41
      2  951 1400 (N31W04 L=259)   791    240        3 ( 892/1284-1008/1495) -90.0 W02-W06 N22-N38
      3 1667 1387 (N12W52 L=307)   396    106        0 (1624/1359-1709/1410) -45.0 W49-W54 N09-N14
      4  426 1105 (N23E38 L=217)  3114    859        8 ( 204/ 954- 720/1252) -32.0 E58-E18 N09-N35
      5 1708  876 (S18W41 L=296)  1944    355        2 (1609/ 732-1759/1053)  43.3 W31-W48 S24-S08
      6  814  937 (N06E14 L=241)   341     75        1 ( 774/ 910- 850/ 970)   0.0 E16-E11 N04-N08
      7  226  816 (N08E62 L=193)   401    143        3 ( 192/ 797- 258/ 832)  -7.1 E65-E57 N07-N09
      8  802  737 (S05E19 L=236)   467    149        1 ( 720/ 703- 875/ 766)   0.0 E24-E14 S06-S04
      9 1192  594 (S22W01 L=255)   249     72        1 (1169/ 552-1209/ 624) -90.0 W00-W01 S25-S19
     10 1541  494 (S40W25 L=279)   360     82        0 (1502/ 470-1578/ 519)  -9.5 W21-W27 S40-S38
     11  554  391 (S23E48 L=206)   441     78        1 ( 518/ 352- 596/ 422) -59.0 E50-E46 S26-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