Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  945 1806 (N40E36 L=347)   571     77        1 ( 906/1784- 982/1828)  18.4 E38-E32 N37-N41
      2  546 1585 (N13E48 L=335)  1906    295        3 ( 500/1464- 587/1732) -58.6 E56-E43 N05-N23
      3  468 1631 (N15E58 L=325)  1096    170        0 ( 442/1552- 508/1711) -54.0 E62-E52 N09-N20
      4 1361 1565 (N36W04 L= 27)   451     72        1 (1340/1533-1383/1604) -38.7 W01-W06 N33-N38
      5  629 1479 (N09E37 L=346)  1681    299        3 ( 528/1389- 750/1564)  61.9 E41-E32 N01-N16
      6 1379 1485 (N31W08 L= 31)   392     60        0 (1347/1467-1410/1505)  11.3 W05-W10 N29-N32
      7 1036 1399 (N15E10 L= 13)   334    136        0 ( 990/1345-1084/1441) -20.6 E13-E05 N13-N16
      8  204 1242 (S14E65 L=318)  2370    241        4 ( 125/1128- 263/1326)  81.3 E83-E60 S20-S07
      9 1325 1145 (N09W13 L= 36)   264     51        1 (1310/1112-1338/1179) -63.4 W11-W14 N07-N11
     10  424 1125 (S16E40 L=343)   264     48        0 ( 400/1104- 446/1149)  45.0 E41-E38 S17-S14
     11  185  978 (S30E63 L=321)  1241    152        0 ( 128/ 927- 237/1026)  24.0 E73-E55 S33-S25
     12  160  919 (S34E67 L=316)   558    104        0 ( 122/ 885- 204/ 948)  14.0 E76-E59 S36-S31
     13  716  774 (S29E11 L= 12)   613    161        0 ( 641/ 748- 786/ 804)  36.9 E15-E07 S31-S25
     14 1406  706 (S12W32 L= 55)   220    108        0 (1355/ 687-1458/ 723)   8.1 W28-W35 S12-S10
     15 1227  576 (S25W27 L= 50)   356    149        0 (1162/ 565-1309/ 588)  24.0 W22-W31 S26-S21
     16 1498  525 (S18W46 L= 70)   959    150        1 (1431/ 502-1572/ 561)  45.0 W41-W50 S21-S14
     17 1112  409 (S38W29 L= 52)  1872    359        2 ( 968/ 334-1257/ 485)  64.7 W23-W33 S47-S28
     18  935  380 (S46W19 L= 42)   320     56        2 ( 908/ 353- 954/ 407) -53.1 W16-W21 S48-S44
  

    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