Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1484 1783 (N54W65 L=162)   774     66        0 (1465/1750-1503/1812)  11.3 W57-W82 N51-N56
      2  605 1682 (N50E40 L= 56)   458     61        2 ( 583/1653- 632/1707)  78.7 E43-E38 N47-N52
      3 1534 1601 (N39W48 L=145)   610    125        0 (1465/1581-1585/1616)  -9.5 W41-W54 N37-N40
      4  899 1437 (N31E08 L= 89)  2439    600        3 ( 650/1289-1109/1554) -28.5 E29-W06 N20-N39
      5 1445 1282 (N18W29 L=126)   504    222        0 (1416/1203-1507/1352)  52.1 W26-W35 N13-N22
      6 1252 1066 (N06W14 L=111)   420    133        1 (1205/1034-1313/1107)  23.2 W11-W18 N04-N08
      7 1056  689 (S17W01 L= 98)   884    209        3 (1001/ 593-1114/ 763)  61.4 E02-W04 S23-S12
      8 1682  524 (S32W54 L=150)  5916    702       11 (1445/ 381-1856/ 739)  31.7 W32-W72 S42-S18
      9  359  514 (S29E58 L= 39)   668    102        1 ( 318/ 479- 400/ 546)  18.4 E66-E51 S31-S26
     10 1309  468 (S34W20 L=116)  1465    253        4 (1216/ 393-1412/ 542)  34.7 W12-W26 S39-S28
     11  823  449 (S33E17 L= 80)   470    103        1 ( 779/ 418- 866/ 485)  29.7 E20-E13 S35-S30
     12 1177  371 (S41W10 L=107)   322     56        1 (1148/ 352-1202/ 393) -31.0 W07-W12 S42-S38
  

    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