Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  614 1806 (N30E56 L= 31)  1643    289        4 ( 582/1740- 649/1873) -50.2 E61-E50 N25-N37
      2 1588 1469 (N36W24 L=111)  2229    375        4 (1451/1325-1708/1598) -23.5 W11-W35 N30-N42
      3  428 1537 (N08E54 L= 34)  1520    238        2 ( 349/1480- 510/1569)   3.6 E61-E45 N05-N10
      4  272 1527 (N05E71 L= 16)   944    197        1 ( 246/1502- 300/1559) -14.0 E78-E66 N03-N07
      5  888 1436 (N13E19 L= 68)   421    136        0 ( 831/1414- 975/1469)  45.0 E22-E14 N10-N16
      6  323 1213 (S13E50 L= 37)  3599    660        4 ( 131/1014- 503/1467)  58.5 E64-E42 S26-N05
      7 1285 1260 (N13W08 L= 95)  1605    403        1 (1143/1066-1384/1391) -28.8 E03-W18 N04-N17
      8  241 1267 (S11E60 L= 28)   729    128        1 ( 199/1209- 278/1319)  60.3 E62-E57 S14-S07
      9  537 1096 (S15E31 L= 57)   382    132        2 ( 489/1036- 596/1154)  71.6 E33-E28 S19-S10
     10 1552  452 (S22W53 L=141)   742    140        0 (1494/ 436-1613/ 470)  21.8 W48-W58 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