Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1141 1584 (N40W05 L=340)   380     89        0 (1092/1566-1181/1600)  15.9 W01-W08 N38-N41
      2  676 1470 (N28E28 L=307)   564    194        0 ( 600/1431- 773/1500)  -7.6 E35-E20 N26-N29
      3 1643 1397 (N29W45 L= 21)  1250    254        1 (1568/1305-1717/1481) -47.7 W40-W51 N23-N34
      4 1430 1424 (N30W27 L=  2)   327     62        1 (1397/1402-1467/1450) -11.3 W24-W29 N28-N31
      5  484 1313 (N16E39 L=296)  1319    282        2 ( 412/1226- 584/1409) -30.5 E48-E30 N11-N21
      6  954 1307 (N19E06 L=329)   688    216        1 ( 864/1247-1029/1359) -28.6 E12-E01 N16-N22
      7  773  726 (S18E14 L=321)   322     92        1 ( 728/ 703- 803/ 747)  21.8 E17-E12 S19-S16
      8 1439  697 (S16W30 L=  5)   320     62        2 (1412/ 667-1469/ 724)  45.0 W27-W31 S17-S13
      9 1701  654 (S18W53 L= 28)   336     46        0 (1684/ 629-1723/ 673)  56.3 W51-W55 S19-S16
     10  291  542 (S35E64 L=271)   705     85        0 ( 248/ 524- 331/ 567)  10.3 E78-E56 S37-S32
  

    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