Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  887 1379 (N30E03 L=  5)  1390    250        3 ( 747/1324-1029/1434)  -5.7 E13-W07 N27-N32
      2 1513 1291 (N14W36 L= 44)   706    155        2 (1434/1258-1575/1315) -32.0 W30-W40 N10-N16
      3 1821 1032 (S08W56 L= 65)  1138    279        0 (1775/ 892-1855/1147)  49.4 W49-W63 S15-S01
      4  835 1105 (N14E10 L=358)   236     75        0 ( 792/1084- 867/1121) -26.6 E13-E08 N12-N15
      5  340 1055 (N16E46 L=322)   774    167        0 ( 291/1018- 397/1115) -33.7 E50-E41 N13-N20
      6  921  831 (S04E09 L=359)  1411    417        1 ( 877/ 688- 958/1018)  81.0 E12-E06 S12-N07
      7 1201  857 (S06W08 L= 16)   360    105        1 (1144/ 823-1240/ 879) -26.6 W04-W10 S08-S04
      8 1887  809 (S24W63 L= 71)   921    258        2 (1846/ 778-1929/ 833) -14.0 W57-W69 S25-S22
      9  863  645 (S15E16 L=352)   225     72        0 ( 829/ 618- 893/ 671)  21.8 E18-E13 S15-S13
     10 1106  540 (S25E02 L=  6)   541    127        0 (1075/ 483-1138/ 605)  54.5 E05-W00 S28-S21
     11  983  464 (S28E12 L=356)   283    114        0 ( 921/ 447-1040/ 494)   0.0 E16-E07 S29-S27
     12  807  321 (S37E30 L=338)  1256    294        2 ( 749/ 253- 853/ 421)  49.4 E37-E25 S42-S28
  

    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