Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1264 1664 (N39W05 L=238)  1450    406        0 (1098/1577-1431/1759) -14.5 E09-W18 N35-N43
      2  360 1625 (N22E65 L=168)  1254    111        4 ( 317/1587- 398/1664)  63.4 E70-E59 N19-N25
      3  497 1126 (S08E34 L=199)   283     81        0 ( 473/1087- 531/1154)  59.0 E35-E32 S10-S05
      4  489  885 (S22E33 L=200)   873    250        1 ( 434/ 812- 537/ 974)  70.0 E37-E29 S27-S16
      5  139  833 (S28E72 L=161)   649     49        1 ( 108/ 820- 170/ 853)   3.8 E81-E66 S28-S26
      6  776  802 (S23E12 L=221)   481    197        0 ( 691/ 750- 860/ 846)  38.7 E17-E06 S27-S19
      7 1267  785 (S16W19 L=252)   318    108        1 (1213/ 760-1300/ 818)  51.3 W15-W20 S18-S13
      8 1275  665 (S23W22 L=256)  3099    861        4 ( 950/ 549-1664/ 818)  33.7 W03-W45 S35-S06
      9 1702  640 (S15W55 L=288)   690     87        0 (1663/ 599-1733/ 679)  71.6 W52-W57 S18-S12
     10  822  587 (S36E06 L=227)   453    126        0 ( 768/ 555- 877/ 628)  -6.3 E10-E01 S37-S34
     11 1727  508 (S22W67 L=300)  1873    130        3 (1668/ 452-1765/ 566) -31.0 W59-W75 S24-S18
     12 1640  461 (S27W60 L=294)   476     85        1 (1618/ 421-1662/ 501) -51.3 W58-W62 S29-S24
  

    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