Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1186 1868 (N51W54 L=262)   453     48        0 (1172/1838-1203/1893)  20.0 W48-W59 N48-N52
      2  592 1609 (N52E16 L=192)   582     99        1 ( 554/1572- 632/1645)   6.3 E20-E11 N49-N54
      3  740 1605 (N47E03 L=206)   364     87        0 ( 703/1588- 782/1631) -38.7 E05-W00 N44-N49
      4  862 1457 (N33W02 L=210)   278     90        1 ( 831/1421- 887/1494)  31.0 E01-W04 N31-N34
      5  958 1195 (N15W00 L=209)   362    157        1 ( 919/1133-1009/1251) -83.7 W00-W01 N10-N19
      6 1878  937 (S24W59 L=268)  1159    106        1 (1850/ 885-1905/ 991)  59.0 W55-W63 S26-S20
      7 1874  859 (S29W59 L=267)   577     79        2 (1839/ 832-1910/ 883)  -4.8 W53-W65 S30-S27
      8 1544  705 (S28W24 L=233)   342    139        4 (1489/ 652-1590/ 746)  20.6 W19-W28 S29-S25
      9 1570  331 (S60W28 L=237)  1255    106        2 (1520/ 308-1616/ 357) -25.3 W20-W39 S64-S55
  

    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