Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  782 1682 (N26E35 L=287)   698     90        1 ( 753/1642- 827/1729) -12.5 E39-E30 N23-N28
      2  542 1683 (N21E53 L=269)   938    108        2 ( 492/1645- 584/1725)  74.1 E55-E49 N17-N24
      3 1198 1376 (N17W01 L=323)   395     58        0 (1176/1339-1225/1409)  90.0 E01-W02 N14-N18
      4 1805 1218 (N27W48 L= 10)  1132    209        1 (1768/1118-1837/1311) -55.0 W43-W52 N21-N32
      5 1080  884 (S13W07 L=329)   201     51        0 (1062/ 859-1110/ 903)  71.6 W05-W08 S15-S12
      6  876  649 (S32W01 L=323)   212     65        0 ( 844/ 626- 909/ 665)  53.1 W00-W03 S34-S30
      7 1112  544 (S32W20 L=342)   580    134        1 (1065/ 508-1164/ 599) -16.7 W14-W24 S33-S29
      8  583  562 (S46E19 L=303)   255     39        0 ( 559/ 544- 606/ 580) -14.0 E21-E16 S47-S45
      9 1521  366 (S28W60 L= 22)   508     71        0 (1485/ 349-1560/ 389)   5.7 W54-W64 S29-S26
     10  718  316 (S61W02 L=323)  1577    334        3 ( 583/ 240- 865/ 385)   0.0 E17-W23 S63-S57
  

    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