Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  723 1730 (N32E43 L=352)  1254    270        3 ( 675/1665- 779/1812) -35.5 E49-E35 N28-N38
      2  842 1746 (N36E36 L=359)   361    105        0 ( 815/1711- 874/1779) -26.6 E39-E31 N34-N38
      3  687 1546 (N19E35 L=359)  1559    393        2 ( 577/1392- 742/1624) -35.2 E45-E28 N09-N23
      4 1141 1541 (N30E07 L= 27)   647    208        0 (1051/1468-1223/1603) -15.9 E14-W00 N27-N33
      5  487 1436 (N08E45 L=350)   437    168        0 ( 462/1382- 535/1490) -29.1 E48-E39 N05-N10
      6 1622 1428 (N37W31 L= 65)   564    129        0 (1569/1378-1680/1461) -11.3 W25-W36 N35-N38
      7 1339 1370 (N25W10 L= 45)   501    142        1 (1277/1327-1421/1408)   0.0 W05-W16 N23-N25
      8  849 1357 (N11E18 L= 16)   283     98        0 ( 824/1310- 878/1404) -33.7 E21-E15 N09-N13
      9 1761 1292 (N32W44 L= 79)   387     79        1 (1730/1269-1785/1310) -11.3 W41-W46 N30-N33
     10  934 1189 (N04E09 L= 25)   291    161        1 ( 873/1150- 964/1266)  76.0 E11-E07 N00-N08
     11 1743 1027 (N15W46 L= 80)   359     89        3 (1716/ 999-1778/1046)  18.4 W43-W49 N13-N16
     12  883  489 (S39W06 L= 40)   290     46        1 ( 850/ 469- 911/ 513)   0.0 W03-W08 S40-S37
  

    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