Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1185 1752 (N49E09 L= 53)  1934    275        6 (1097/1666-1253/1828) -15.5 E18-W00 N42-N56
      2  402 1540 (N15E57 L=  6)  1887    254        5 ( 367/1429- 448/1646) -47.1 E67-E52 N08-N22
      3  628 1198 (S00E27 L= 35)  1147    297        2 ( 545/1081- 713/1310)  68.2 E30-E24 S08-N07
      4 1385 1262 (N19W17 L= 79)   213     98        0 (1338/1237-1426/1296)  -9.5 W13-W19 N18-N20
      5  957 1095 (N00E05 L= 57)   183    109        0 ( 943/1040- 965/1144) -71.6 E06-E04 S03-N03
      6 1208 1005 (N00W11 L= 73)   206    131        0 (1197/ 939-1228/1063) -74.1 W10-W12 S03-N04
      7  852  989 (S08E09 L= 53)   241    104        0 ( 807/ 951- 895/1019)  51.3 E11-E07 S10-S05
      8 1799  770 (N00W60 L=122)  3298    266        7 (1738/ 658-1855/ 903) -27.8 W50-W72 S04-N07
      9 1462  848 (S03W30 L= 92)   517    105        2 (1418/ 810-1504/ 888)  51.3 W27-W32 S05-N00
  

    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