Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  899 1786 (N61W15 L=288)   623    100        1 ( 849/1762- 947/1808)   0.0 W08-W22 N58-N63
      2 1384 1715 (N40W51 L=323)   496     79        0 (1347/1686-1418/1753)   8.7 W45-W58 N38-N41
      3  445 1624 (N58E48 L=224)   477     67        0 ( 416/1605- 474/1643)  14.0 E51-E43 N56-N59
      4 1068 1590 (N41W18 L=290)   683    161        0 (1033/1535-1141/1625)  10.3 W13-W24 N37-N43
      5 1237 1585 (N36W30 L=302)   259     50        0 (1208/1568-1269/1600) -26.6 W27-W32 N34-N36
      6 1732 1520 (N18W68 L=341)  1184    108        0 (1691/1485-1786/1545)   0.0 W61-W81 N16-N19
      7 1619  690 (S27W33 L=306)   552     96        1 (1599/ 642-1635/ 740)  63.4 W31-W34 S29-S23
      8 1431  683 (S22W19 L=291)   297     99        0 (1408/ 642-1458/ 728)  45.0 W16-W21 S24-S19
  

    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