Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  770 1541 (N23E23 L=354)   369    154        0 ( 727/1494- 797/1616)  82.9 E25-E22 N20-N29
      2  592 1443 (N15E33 L=345)   936    153        5 ( 545/1375- 642/1519) -48.4 E37-E28 N11-N20
      3 1427 1423 (N24W21 L= 39)   287    119        0 (1390/1373-1461/1472) -45.0 W18-W23 N21-N26
      4  968 1418 (N18E08 L=  9)   256     72        0 ( 928/1393- 997/1439) -26.6 E10-E06 N16-N18
      5  605 1309 (N07E29 L=348)   408     70        1 ( 574/1277- 628/1347)  90.0 E31-E27 N05-N09
      6  538 1246 (N03E33 L=345)   507     87        2 ( 492/1212- 578/1277)  21.8 E35-E30 N00-N04
      7  182  853 (S23E64 L=313)   565     53        1 ( 156/ 824- 203/ 883)  38.7 E67-E60 S24-S20
      8  294  744 (S29E53 L=325)   808    261        0 ( 263/ 632- 316/ 869)  56.3 E60-E48 S36-S21
      9 1736  506 (S25W67 L= 85)  1413    153        2 (1665/ 473-1805/ 534)  18.4 W59-W79 S28-S22
     10 1319  496 (S34W29 L= 47)   420    103        0 (1273/ 469-1367/ 528)  45.0 W26-W31 S36-S31
     11  947  321 (S53W06 L= 23)   580     39        1 ( 917/ 292- 978/ 350)  31.0 W02-W08 S55-S50
  

    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