Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1179 1600 (N35W29 L=156)  1901    437        3 (1010/1545-1385/1637) -17.8 W16-W44 N28-N39
      2  963 1472 (N34W09 L=136)   503    116        0 ( 945/1422- 982/1525)  50.2 W06-W12 N30-N37
      3  967 1373 (N28W06 L=133)   253     44        0 ( 938/1358-1000/1390) -36.9 W04-W08 N26-N29
      4  330 1017 (N22E47 L= 80)   334     90        0 ( 297/ 987- 368/1050) -39.8 E49-E43 N19-N24
      5 1860  952 (S22W57 L=183)   945    110        1 (1821/ 925-1900/ 981) -31.0 W52-W63 S24-S18
      6 1538  796 (S20W25 L=151)   282     70        0 (1516/ 768-1557/ 842)  36.9 W22-W27 S21-S17
      7 1644  715 (S29W32 L=158)   936    212        1 (1560/ 670-1736/ 777)  -3.4 W24-W41 S30-S26
      8 1526  579 (S34W20 L=147)   861    131        2 (1485/ 533-1563/ 642)  39.8 W16-W24 S36-S30
      9  458  411 (S17E58 L= 68)   495    110        1 ( 431/ 367- 494/ 451)  81.9 E60-E55 S20-S13
     10  401  383 (S18E68 L= 58)  1275    136        2 ( 354/ 353- 432/ 415) -45.0 E75-E63 S20-S15
  

    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