Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  958 1728 (N37E28 L=299)  1137    220        3 ( 916/1625-1016/1809) -31.4 E36-E18 N31-N42
      2  661 1682 (N26E45 L=282)   889    166        2 ( 619/1613- 708/1751) -23.6 E54-E38 N22-N29
      3 1518 1699 (N57W19 L=347)  1614    190        1 (1439/1654-1586/1730)  17.5 W10-W30 N50-N61
      4  185 1000 (S25E63 L=265)  1224    191        2 ( 143/ 915- 245/1061)  39.8 E70-E54 S30-S20
      5  533  747 (S33E25 L=303)   696    204        1 ( 441/ 719- 641/ 767)  23.6 E32-E16 S35-S28
      6 1050  597 (S28W14 L=341)   609    137        0 (1020/ 529-1084/ 657) -71.6 W12-W15 S32-S23
      7 1419  532 (S21W40 L=  7)   317     95        1 (1390/ 493-1451/ 577) -71.6 W38-W41 S23-S17
      8 1268  453 (S30W34 L=  2)   356    136        0 (1229/ 421-1321/ 507) -10.3 W28-W39 S31-S27
      9  477  398 (S62E32 L=295)  1876    258        5 ( 394/ 336- 559/ 472) -19.3 E48-E20 S65-S56
     10  691  330 (S59E02 L=325)   343     47        1 ( 668/ 309- 713/ 351)   0.0 E04-W01 S61-S57
     11  892  214 (S61W28 L=356)   603     59        0 ( 869/ 180- 921/ 245) -15.9 W22-W36 S63-S58
     12 1231  201 (S47W56 L= 24)   809     97        1 (1196/ 165-1275/ 245)  -5.7 W47-W67 S48-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