Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1001 1792 (N52W28 L= 20)   553     55        1 ( 962/1774-1037/1813) -38.7 W24-W30 N49-N54
      2 1568 1685 (N27W63 L= 55)  1092    111        2 (1516/1655-1609/1711) -15.3 W58-W70 N24-N29
      3 1501 1593 (N24W49 L= 41)   706    166        0 (1414/1571-1572/1611) -31.0 W43-W54 N20-N26
      4 1010 1509 (N32W12 L=  5)   746    202        0 ( 913/1482-1109/1532) -26.6 W05-W17 N27-N35
      5 1699 1481 (N12W59 L= 51)  1514    159        4 (1635/1414-1745/1532) -77.5 W54-W63 N07-N16
      6 1359 1499 (N22W34 L= 26)   545    148        1 (1281/1482-1430/1528) -21.8 W28-W38 N19-N23
      7  430 1068 (N19E36 L=316)   377     74        2 ( 401/1032- 455/1095)  45.0 E38-E34 N17-N20
      8 1705  643 (S37W40 L= 32)   359     44        2 (1680/ 621-1728/ 670)  71.6 W37-W42 S39-S35
      9 1426  502 (S38W13 L=  6)  1061    295        2 (1323/ 407-1530/ 591)  22.4 W04-W22 S43-S34
     10 1520  299 (S61W23 L= 16)   749     86        0 (1476/ 274-1568/ 325)  -8.1 W17-W31 S63-S59
     11  595  277 (S32E57 L=295)  1039     91        1 ( 562/ 248- 635/ 319) -90.0 E61-E54 S34-S28
  

    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