Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1121 1635 (N37W19 L=239)   357     49        0 (1094/1609-1155/1660)   9.5 W16-W22 N34-N38
      2 1509 1635 (N29W49 L=269)   415     67        1 (1473/1606-1536/1656) -63.4 W45-W50 N27-N31
      3  476 1337 (N28E32 L=188)  3877    699        7 ( 176/1241- 761/1429)  -4.9 E68-E10 N22-N35
      4 1682  875 (S19W41 L=261)  2197    607        4 (1523/ 678-1827/1055)  31.3 W27-W55 S28-S09
      5 1278  812 (S16W12 L=232)   317    100        0 (1245/ 770-1322/ 838)  21.8 W09-W14 S17-S14
      6 1469  737 (S24W24 L=244)   381     88        0 (1450/ 690-1489/ 788) -80.5 W22-W25 S26-S20
      7 1260  684 (S23W09 L=229)   320    108        0 (1217/ 647-1310/ 705)   8.1 W05-W12 S25-S22
      8  468  647 (S13E43 L=177)   282     52        1 ( 443/ 624- 499/ 663)  14.0 E44-E39 S13-S11
      9  385  488 (S21E58 L=162)  1077    224        1 ( 312/ 441- 459/ 538)   7.8 E70-E48 S23-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