Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  976 1610 (N39E11 L=126)  1009    169        2 ( 904/1558-1056/1653) -17.1 E16-E03 N36-N41
      2  492 1511 (N27E45 L= 91)  1705    245        5 ( 419/1417- 546/1592) -34.5 E56-E38 N20-N31
      3 1044 1380 (N24E02 L=135)   349     79        0 (1026/1331-1058/1422) -80.5 E03-E01 N20-N26
      4 1810 1222 (N20W58 L=195)   369    100        0 (1792/1176-1828/1266) -68.2 W56-W59 N17-N22
      5  639  777 (S17E23 L=114)  2914    633        2 ( 411/ 591- 814/ 987)  43.1 E43-E12 S32-S02
      6 1302  813 (S09W19 L=156)   338    130        1 (1248/ 777-1336/ 846)  45.0 W16-W21 S11-S06
      7 1631  685 (S14W46 L=182)   407    105        0 (1605/ 661-1663/ 705) -11.3 W43-W48 S15-S12
      8 1604  526 (S25W49 L=186)  1746    330        2 (1481/ 433-1705/ 634)  49.4 W41-W55 S32-S16
      9  292  512 (S40E66 L= 71)  1018     91        0 ( 269/ 472- 325/ 546) -90.0 E85-W00 S44-S37
     10 1400  424 (S34W36 L=173)   438     70        1 (1364/ 409-1437/ 442)   9.5 W32-W39 S34-S32
  

    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