Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  853 1730 (N31E34 L=139)  1601    330        0 ( 779/1621- 947/1787)  90.0 E38-E26 N21-N36
      2 1264 1675 (N40E04 L=169)   720    174        0 (1196/1614-1313/1728) -24.4 E10-W01 N36-N42
      3  446 1665 (N18E61 L=111)   483     53        1 ( 418/1649- 475/1684)  36.9 E64-E58 N16-N19
      4  449 1481 (N06E49 L=124)   986    153        2 ( 383/1450- 508/1526)  40.6 E54-E44 N03-N09
      5 1778 1422 (N41W46 L=219)   799    112        1 (1740/1373-1812/1468) -29.7 W41-W50 N37-N43
      6  314  941 (S28E47 L=126)   751    249        1 ( 270/ 828- 354/1052)  70.3 E51-E43 S35-S21
      7  797  897 (S20E10 L=163)   581    194        0 ( 718/ 831- 867/ 939) -10.3 E15-E04 S22-S18
      8  440  781 (S35E34 L=139)   842    257        0 ( 403/ 689- 497/ 872) -63.4 E37-E27 S40-S29
      9 1563  736 (S08W40 L=213)   381     60        0 (1544/ 701-1582/ 770) -78.7 W38-W42 S10-S05
     10 1125  673 (S24W16 L=189)   221     50        2 (1107/ 645-1146/ 699) -56.3 W14-W17 S25-S22
     11 1404  461 (S27W43 L=216)   991    267        0 (1280/ 431-1516/ 500)  31.0 W34-W49 S30-S21
     12 1098  358 (S44W30 L=203)   312     63        0 (1073/ 334-1125/ 389)  90.0 W28-W31 S45-S40
     13  652  323 (S64E07 L=166)  1321    150        3 ( 581/ 275- 706/ 370)  -6.3 E17-W01 S67-S60
     14 1558  341 (S28W67 L=239)  1656    157        2 (1474/ 307-1628/ 364)  13.4 W57-W79 S31-S25
     15 1332  297 (S39W52 L=225)   470     45        1 (1308/ 274-1359/ 325) -15.9 W48-W55 S40-S37
     16 1160  294 (S46W40 L=213)   328     32        1 (1143/ 268-1180/ 317) -31.0 W37-W42 S47-S44
     17  755  292 (S62W07 L=180)   331     21        0 ( 735/ 272- 774/ 312)  76.0 W05-W09 S64-S60
     18 1263  256 (S44W52 L=224)  1721    193        4 (1194/ 221-1347/ 304)  10.0 W41-W60 S46-S40
  

    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