Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  370 1567 (N29E63 L=212)  2148    294        3 ( 333/1459- 409/1665) -24.1 E82-E52 N22-N35
      2  562 1428 (N22E37 L=238)   350     52        1 ( 542/1395- 582/1458)  76.0 E38-E35 N19-N24
      3 1053  981 (S01W02 L=277)   546    129        1 ( 994/ 956-1101/1035)   9.5 E01-W05 S02-N02
      4 1553  969 (N03W35 L=310)   591    104        0 (1523/ 920-1583/1019)  59.0 W32-W36 N00-N05
      5  178  844 (S18E67 L=208)  2022    156        4 ( 153/ 770- 207/ 923)  55.0 E73-E61 S23-S13
      6  925  825 (S12E04 L=271)   307     66        1 ( 883/ 799- 960/ 841)   0.0 E07-E02 S12-S10
      7 1734  695 (S13W56 L=331)  2238    316        5 (1682/ 605-1811/ 802)  72.9 W51-W62 S18-S05
      8  826  757 (S17E10 L=265)   231     79        1 ( 783/ 735- 847/ 783)  26.6 E13-E08 S18-S14
      9 1103  717 (S17W08 L=283)   436    166        1 (1028/ 694-1162/ 769) -12.5 W02-W11 S18-S14
     10  336  710 (S25E49 L=226)   450    107        1 ( 289/ 692- 379/ 737)  21.8 E54-E44 S26-S22
     11 1502  664 (S17W36 L=311)   790    248        0 (1447/ 623-1558/ 734)  49.4 W31-W39 S19-S11
     12  951  497 (S33W00 L=275)   386     66        0 ( 918/ 478- 985/ 517)  26.6 E02-W02 S34-S31
  

    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