Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  676 1578 (N35E22 L=325)   272     47        1 ( 653/1557- 701/1602)  56.3 E23-E20 N33-N36
      2  657 1413 (N23E22 L=325)   240     74        0 ( 641/1382- 681/1451) -59.0 E23-E20 N21-N26
      3 1699 1351 (N15W48 L= 35)   590    139        1 (1652/1290-1742/1407) -66.8 W44-W51 N11-N18
      4 1480 1277 (N11W30 L= 17)   367     95        1 (1458/1241-1507/1328) -71.6 W28-W31 N08-N14
      5 1052 1174 (N06W02 L=349)   989    254        1 ( 944/1064-1128/1255) -50.7 E03-W06 N00-N11
      6  638 1123 (N06E23 L=324)   190     56        0 ( 607/1103- 670/1143) -26.6 E25-E21 N04-N06
      7  648 1066 (N02E23 L=324)   473    154        4 ( 603/1036- 704/1097) -26.6 E25-E19 N00-N04
      8 1665  842 (S16W42 L= 29)   596    119        2 (1631/ 781-1694/ 896)  66.8 W39-W44 S19-S12
      9  198  860 (S07E61 L=286)   704    151        2 ( 165/ 815- 238/ 894)  24.0 E65-E55 S09-S05
     10  898  645 (S25E10 L=337)   203     60        1 ( 867/ 618- 931/ 669) -53.1 E12-E08 S26-S22
     11 1443  448 (S41W30 L= 17)  1461    346        2 (1291/ 373-1604/ 521)  16.3 W18-W42 S46-S36
  

    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