Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1328 1571 (N32W21 L=  2)  1668    299        3 (1219/1466-1443/1690) -50.9 W14-W28 N24-N40
      2  515 1525 (N28E37 L=305)   318     27        0 ( 494/1507- 541/1541)  14.0 E38-E34 N26-N29
      3  559 1371 (N18E31 L=311)   529     99        0 ( 537/1318- 576/1422) -63.4 E32-E29 N14-N20
      4  424 1245 (N10E39 L=302)   673    170        2 ( 355/1225- 495/1271)  -5.2 E45-E34 N08-N11
      5  750 1214 (N08E17 L=324)   298    107        0 ( 715/1159- 793/1269)  60.3 E18-E14 N04-N11
      6 1082 1080 (S00W03 L=345)   290     49        0 (1048/1058-1113/1101) -14.0 W01-W05 S01-N01
      7 1039  898 (S11W01 L=342)   649    197        0 (1000/ 802-1073/ 999)  74.7 E01-W02 S16-S05
      8 1855  801 (S14W62 L= 43)  1425    243        5 (1800/ 758-1900/ 850)  20.0 W56-W67 S17-S11
      9  731  658 (S26E19 L=322)   546    129        0 ( 698/ 599- 759/ 716)  60.3 E22-E17 S29-S22
     10  414  650 (S26E43 L=298)   849    145        0 ( 367/ 606- 446/ 715)  41.2 E48-E40 S28-S21
     11 1011  546 (S33E00 L=341)   684    165        1 ( 946/ 509-1076/ 601)  36.9 E05-W03 S35-S29
  

    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