Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1196 1747 (N44E15 L=219)  1085    160        2 (1149/1674-1251/1799) -77.0 E16-E12 N37-N50
      2 1278 1763 (N49E10 L=224)   508     50        1 (1258/1732-1300/1789) -38.7 E12-E07 N46-N52
      3 1554 1432 (N34W22 L=255)   990    188        3 (1526/1346-1582/1520) -85.6 W18-W24 N28-N41
      4 1301 1432 (N25W04 L=238)   219     64        0 (1272/1414-1337/1453)   0.0 W01-W06 N24-N26
      5 1740 1348 (N36W40 L=274)  1604    300        1 (1666/1209-1797/1448) -36.3 W32-W47 N29-N41
      6  764  780 (S27E08 L=225)   378    146        1 ( 709/ 734- 827/ 819)  49.4 E10-E04 S30-S22
      7 1201  704 (S18W20 L=253)   695    268        2 (1116/ 666-1300/ 775)  18.4 W13-W25 S20-S14
      8  527  652 (S42E23 L=211)   254     56        0 ( 506/ 625- 544/ 678) -56.3 E25-E21 S43-S39
      9  983  485 (S37W16 L=250)   522    126        1 ( 942/ 446-1038/ 525)  76.0 W13-W17 S41-S33
     10 1127  352 (S41W34 L=267)   309     67        0 (1099/ 335-1155/ 374)  76.0 W31-W35 S42-S38
     11 1319  233 (S41W60 L=294)   405     56        0 (1293/ 214-1349/ 259)   0.0 W54-W65 S41-S39
  

    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