Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1744 1515 (N32W61 L=106)   815    146        2 (1708/1486-1776/1539)   0.0 W56-W64 N30-N33
      2 1568 1516 (N34W44 L= 89)   494     57        1 (1539/1493-1598/1538) -36.9 W41-W46 N32-N35
      3  395 1488 (N33E50 L=355)   983    251        2 ( 338/1436- 456/1534) -23.2 E57-E42 N30-N36
      4  733 1386 (N27E20 L= 25)   267     75        2 ( 697/1370- 778/1401)   0.0 E22-E16 N25-N27
      5 1742 1290 (N18W53 L= 98)   717    208        2 (1667/1212-1801/1362) -45.0 W47-W58 N12-N22
      6  257 1248 (N17E57 L=348)   849    151        2 ( 238/1191- 283/1302) -66.8 E59-E54 N13-N20
      7 1645 1211 (N13W43 L= 88)   677    162        2 (1599/1127-1680/1289) -73.3 W40-W45 N08-N18
      8  912 1188 (N14E07 L= 38)   299    109        0 ( 869/1144- 960/1240) -45.0 E09-E03 N11-N17
      9 1319  970 (S00W18 L= 63)   196     73        0 (1283/ 951-1358/ 990)   0.0 W15-W20 S01-N01
     10  232  662 (S20E63 L=342)   860    208        1 ( 215/ 605- 251/ 708)  48.8 E67-E60 S24-S16
     11 1209  571 (S25W12 L= 57)   334    100        1 (1160/ 558-1260/ 584)  -8.1 W08-W15 S26-S24
     12  419  546 (S27E47 L=358)   715    148        1 ( 382/ 519- 466/ 569)  21.8 E51-E41 S29-S25
  

    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