Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1612 1635 (N21W61 L= 50)   744    121        0 (1558/1604-1664/1665) -49.4 W58-W64 N17-N24
      2 1430 1617 (N27W46 L= 35)   316     62        1 (1409/1584-1449/1646)  76.0 W43-W47 N24-N28
      3 1309 1573 (N28W36 L= 24)   599    115        1 (1242/1548-1367/1595) -35.5 W31-W38 N25-N30
      4  403 1110 (N27E37 L=311)   450    142        1 ( 366/1059- 460/1142)  12.5 E41-E32 N24-N28
      5  410  997 (N20E39 L=310)   383     88        0 ( 369/ 970- 445/1017) -38.7 E41-E36 N17-N21
      6 1817  893 (S25W48 L= 37)  1281    230        1 (1744/ 835-1864/ 983)  17.1 W40-W54 S27-S21
      7 1674  878 (S21W36 L= 24)   800    171        0 (1650/ 820-1725/ 945) -77.5 W33-W38 S25-S16
      8 1159  701 (S15E01 L=347)   436    106        3 (1131/ 638-1199/ 740)  38.7 E04-W01 S18-S12
      9 1339  574 (S28W06 L=355)   841    244        0 (1243/ 489-1424/ 691)  28.1 E02-W14 S31-S23
  

    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