Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  795 1609 (N25E32 L=284)  1529    484        2 ( 682/1429- 896/1801) -28.5 E53-E18 N15-N35
      2 1095 1673 (N38E16 L=300)  1301    222        3 (1035/1575-1152/1758) -42.7 E23-E08 N31-N43
      3 1578 1666 (N57W28 L=343)  1240    130        1 (1522/1631-1639/1691)  29.5 W20-W43 N50-N63
      4  275  945 (S26E51 L=265)  1535    301        1 ( 206/ 827- 382/1019)  35.8 E60-E40 S35-S20
      5  678  687 (S32E13 L=303)   762    224        2 ( 574/ 670- 798/ 705)  21.8 E20-E04 S35-S28
      6 1198  522 (S28W26 L=342)   825    149        2 (1153/ 458-1231/ 588) -71.6 W23-W28 S32-S23
      7 1478  487 (S22W47 L=  2)   278     77        2 (1447/ 456-1508/ 515)  90.0 W44-W48 S24-S19
      8  534  375 (S61E24 L=292)  2396    356        3 ( 403/ 293- 665/ 461)  -6.2 E43-E04 S65-S56
      9 1331  402 (S31W42 L=357)   740    277        0 (1200/ 368-1436/ 441)  14.0 W31-W51 S34-S28
     10  832  276 (S58W15 L=331)   413     64        0 ( 802/ 248- 863/ 302) -12.5 W11-W20 S59-S56
     11 1268  176 (S48W65 L= 21)   620     47        1 (1250/ 148-1282/ 207) -10.0 W57-W74 S49-S45
  

    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