Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1459 1671 (N50W10 L=334)   773    104        1 (1427/1625-1492/1709) -35.5 W06-W14 N46-N53
      2  713 1505 (N14E33 L=292)   705    231        2 ( 651/1431- 753/1609) -81.3 E35-E30 N08-N21
      3 1073 1545 (N26E13 L=312)   274    151        0 (1020/1488-1114/1599) -24.0 E17-E08 N24-N28
      4 1714 1448 (N43W39 L=  3)  1924    276        3 (1602/1361-1802/1517)  -7.1 W26-W50 N39-N45
      5 1254 1333 (N19W05 L=329)   622    235        1 (1231/1246-1276/1455) -77.0 W01-W07 N13-N26
      6 1560 1190 (N20W27 L=352)   547    247        0 (1486/1078-1598/1258) -30.3 W20-W32 N14-N22
      7 1706 1058 (N17W42 L=  6)   648    135        0 (1686/1001-1721/1130) -60.3 W39-W43 N13-N21
      8 1125  626 (S24W18 L=342)  1249    426        1 (1068/ 534-1251/ 793) -30.5 W11-W29 S28-S14
      9  709  647 (S35E09 L=315)   663    284        0 ( 590/ 633- 854/ 661)  27.9 E17-W00 S39-S30
     10  796  498 (S43W02 L=326)  1129    289        3 ( 681/ 472- 935/ 530)  30.5 E07-W10 S47-S35
     11 1624  435 (S19W64 L= 28)   698     92        2 (1579/ 407-1666/ 461)  51.3 W60-W67 S21-S16
  

    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